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As  we  mentioned  in  the  previous  progress  report,  in  addition  to  the  studies 
on  swirling  flows  of  Newtonian  and  non-Newtonian  fluids,  we  also  studied  other 
closely  related  problems  involving  flows  of  non-Newtonian  fluids  which  can  shed 
light  on  the  mechanics  of  such  fluids.  In  all,  we  studied  the  following  problems: 

(i)  Flow  of  non-Newtonian  fluids  of  the  rate  type  between  two  parallel 
plates  rotating  with  different  angular  speeds  about  a  common  axis. 

(ii)  Flow  of  non-Newtonian  fluids  of  the  integral  type  between  parallel 
plates  rotating  about  distinct  axes. 

(iii)  Flow  of  non-Newtonian  fluids  due  to  torsional  and  longitudinal 
oscillations. 

(iv)  Flow  of  non-Newtonian  fluids  in  pipes  of  varying  cross-sections. 

(v)  Flow  of  shear  thinning  fluid  between  intersecting  planes. 

The  results  of  our  work  have  been  published  in  archival  journals  and  copies 
of  these  have  been  attached  to  this  final  report.  Here,  we  discuss  briefly  our 
contributions  to  the  understanding  of  each  of  the  above  problems.  It  would  be  fair 
to  say  that  this  study  has  contributed  significantly  to  the  understanding  of 
swirling  flows  in  both  Newtonian  and  non-Newtonian  fluids  and  this  is  reflected 
in  the  various  invitations  we  have  received  both  nationally  and  internationally  to 
present  our  work. 

The  grant  also  supported  the  doctoral  dissertation  of  Mr.  Z.H.  Ji  which  was 
titled  "Multiplicity  of  von  Kami  an  flows  of  viscoelastic  fluids"  awarded  in  1990. 
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atical  preliminaries 


Let  QqcE  denote  the  reference  configuration  of  the  body.  We  shall  denote 
by  X  €  Q0,  particles  belonging  to  the  body  in  the  reference  configuration.  Let  the 


configuration  of  the  body  at  time  t  be  denoted  by  fl*  and  let  x  denote  the  position  of 

X  at  time  t.  By  the  motion  of  the  body,  we  mean  a  one-to-one  invertible  mapping  x 
through 

*  =  (1) 
The  deformation  gradient  F  associated  with  the  motion  is  defined  through 


F  = 


ft 

dX  ' 


(2) 


We  shall  assume  that  F  is  invertible.  Let  4  denote  the  position  of  the  particle  at 
time  x,  0  <,  x  <,  t.  Then 


$  =  Z(x,  r)  =  z{z\x,t),t)  =  Xl(x,t) .  (3) 

Zt  is  called  the  relative  motion  The  relative  deformation  gradient  Ft  (x)  is  defined 

through 


*)(*)  = 


K 

dx 


The  history  of  the  stretch  tensor  Ct  (x)  is  defined  through 


(4) 


C,(T)  =  /f(rtf(T). 

The  velocity  of  the  particle  v(x,  t)  is  defined  through 


dx 


(5) 


(6) 


and  the  velocity  gradient  L  is  defined 


dv 

L(x,t)  =  grad  v(x,t)  =  — . 


(7) 


It  immediately  follows  that 

L(x,t)  =  FF~l.  (8) 

We  next  record  the  basic  balance  laws. 

The  conservation  of  mass  takes  the  form 

(9) 

where  p  denotes  the  density  of  the  fluid.  Since  we  shall  be  concerned  with 
incompressible  materials,  we  shall  require  the  constraint 


& 

< 

II 

o 

•* 

(10) 

or 

det  F  =  1 . 

(11) 

The  balance  of  linear  momentum  is  given  by 

div  T +pb^p~, 
dt 


(12) 


where  T  denotes  the  Cauchy  stress,  b  the  external  body  force  field  and  d/dt  the 
usual  material  time  derivative. 


RESEARCH  CARRIED  PITT 

(i)  Ejow_of  fluids  of  the  rate  type  between  two  parallel  nlates  rotating  with 
differing  angular  speeds  ahmit  a  mmmnn  aria 


Amongst  the  many  models  for  the  fluids  of  the  rate  type,  one  that  is  very 
popular  amongst  rheologists,  and  includes  the  classical  Navier-Stokes  model  and 


i 


the  Maxwell  model  as  special  cases,  is  the  Oldroyd-B  fluid. 

The  constitutive  equation  for  the  extra  stress  S  for  the  Oldroyd-B  fluid  is 
given  by  (cf.  Oldroyd  [1]) 

S  =  pl  +  7\ 

S  +  A ,(S  -  IS  -  SLt)  =  44,  +  A44  -  LA,  -  44)]  (13) 

where  the  overscore  dot  denotes  the  material  time  derivative,  p  is  the  viscosity, 
and  At  and  A2  are  material  constants,  referred  to  as  the  relaxation  time  and  the 

retardation  time,  respectively.  Also, 

L  =  grad  v, 

\  =  L  +  lJ, 

and 

(  )=  jj ;+[grad  (  )]v. 

We  are  concerned  here  with  the  flow  that  is  induced  by  the  rotation  of  two 
infinite  disks.  The  disks  are  parallel  to  one  another  and  are  located  at  z  =  0  and 
z  «  /  of  a  cylindrical  polar  coordinate  system.  The  disks  rotate  with  angular 
velocity  and  Q2,  respectively.  The  separation  between  the  disks  is  d  and  the 

fluid  velocity  v  has  radial,  azimuthal  and  axial  components  designated  by  u,  v, 
and  w. 

We  assume  a  velocity  field  of  the  form  (cf.  von  Karman  [2]): 

{u,v,w}  =  [rF',rG,- IF }  (17) 


(14) 

(15) 

(16) 


where  F  and  G  are  functions  of  the  axial  coordinate  z  only.  A  lengthy  but  straight 
forward  manipulation  yields  the  following  equations  (cf.  Ji,  Rajagopal  &  Szeri  [3]): 


5,  +  A,[r(F'5^ 


S^  +  A^F'S*., 


Sn  + \[r(F'S„.r 


Se0  +  A,[r(F  Sgg  r 


Set  +  Aj[r(F  Sg,  r 


2F"Sj)  -  2F’S„  -  2FSni  +  fl] 

:  /x{2  F'  -  2A2[r2F"2  -  2(F'2+FF")]| 


-  2  F"Sie  -  G'S„)  -  2  F'S*  -  2FS„.,  +  GS„B] 
=  -2/i  A2r2F"G' 

2F"StI)  +  F'Sn  -  2FSnt  +  G5raS] 

=  /ir[F"  +  2  A2  (3F'F"  -  FF~)] 


-  2G'S.)  -  2F'S„  -  2FSe,  +  GS„$] 

=  n  {2F'  -  2Aj[r2G'2  +  2(F'2  +  FFW)]} 

■G'S.)  +  F/$h-2FV.  +  GS*#] 

=  a^[G'  +  2  A2(3  F'G'  -  FG")] 


5a  +  Aj{rF'Sn  r+  4F/Sb  -  2FS„,  +  GS„.,} 

=  7?o[-  4F'  +  8A2(FF"  -  2F/2)] 


(18) 


The  stress  components  have  to  satisfy  the  following  set  of  ordinary 
differential  equations 

A  -  2A,(FA'  +  F'A)  =  /i[2F'  -  4A2  (F/2  +  FF")] 

B  -  2A1(F5'  +  F'B)  =  /i[2F'  -  4A2(F'2  +  FF")] 

Q  -  2A,(FQ'  +  G'P)  =  -2/iA2G'2 
F  -  A,(2FF'  -  2F'P  +  G7?)  =  /i[G'  +  2  A2(3F'G'  -  FG")] 

-2Aj(2F/?'  -  2F'R)  =  ^i[-AF'  +  gA^FF'  -  2F'2)] 

-2A,(2 FX'  -  F“Z)  =  -2/i  A2F"2 
K  -  AftFY'  +  G'Z  +  F"P)  =  -2/1  A2F"G' 

Z  -  A1(2FZ/  -  2F'Z  +  F'7?)  =  /i[F'  +  2A2(3  F'F*  -  FFW)] .  (19) 

Here  A,  B,  ...  etc.,  are  defined  through 


(20) 


*>  = 


fA  0 
0  B  0 
,0  0  F, 
r0  0  Z' 
OOF, 


U  P  oj 


X  Y  (T 


s2=  r  q  o  . 


0  0  0; 


form 


The  balance  of  linear  momentum,  in  the  absence  of  body  force,  takes  the 


pr(F'2  -  2FF"  -  G2)  =  ^  ^ 

or  or  dz  r 

2pr(F’G  -  FG’)  =  ^  +  +  7$. 

,  cfc)  oS  S 

4pFF'  =  -- ^  +  — »  +  — “  +  -a . 
oz  or  oz  r 


PxPM+^-r2 

then,  substituting  for  S^., ...,  Szz  we  obtain  the  equations  of  motion  as 


3  X-Q  +  Z'  =  p(F'2  -  2FF”  -G2+K) 

4Y  +  p'  =  2p(F'G-FG') 

p  =  —  r2  +  R{z)  +  2  J  Z(z)dz  -  2pF'2  + constant 


AsB. 

To  non-dimt  nsionalize  the  problem  we  employ  the  transformation: 


z  =  dz,  F  =  Q/(z),  G=Qlg(z),  E=vfQxd2, 

w=ax a„  a2  =  q,a2,  p  =  a2/a1,  k  =  K/nl, 
Q  =  ^Q,  p  =  ^p,  R=pnxR, 

d  d 

X=^X,  Y  =  ^Y,  Z  =  ^Z. 
d2  d2  d 


(25) 


Here  E  is  the  Ekman  number,  W  is  the  Weissenberg  number  and  ^  is  a  measure 
of  retardation  time  relative  to  relaxation  time. 

The  above  set  of  ordinary  differential  equations  read,  alter  dropping  the 
overscore  bars 

Q-lWifQ'  +  g^^-lpWg'1 

P  -  W(2fP'  -  2  fP  +  g'R)  =  g'  +  2pW{SfY  ~  fg") 

R  -  2W(fR'  -2fR)  =  -4/'  +  8 -  2  f'2) 

X  -  2  W(f'  +  f‘ 2)  =  -2f3W  f2"2 

Y  -  W{2fY  +  g'Z  +  fT)  =  -2 pWfY 

Z  -  W(  2  fZ'  -2  f'Z  +  f"R)  =  f"  +  2pw(3fr  -  ffm) 

3x-e+r=-|(/r'-^+f-f) 

4  r -/*%  (26) 

We  next  record  the  boundary  conditions  appropriate  for  the  flow  under 
consideration.  Since  the  fluid  adheres  to  the  boundary 

/( 0)  =  0,  /'(0)  =  0,  *(0)  =  1, 

/(l)  =  0,  /'(1)  =  0,  g(l)  =  s  (27) 

where  s  =  Q2/&i-  As  we  have  eight  differential  equations  governing  the  motion  of 

the  fluid  and  the  above  conditions  yield  but  six  of  the  required  twelve  boundary 
conditions,  we  have  to  augment  (16);  to  this  end  we  evaluate  the  stress 
components  Z,  P  and  R  at  both  z  =  0  and  z  =  1,  as  they  are  determined  in  effect  by 
(16).  Thus 


z(o)=no),  z(i)=ni), 

m=g\o),  p(\)=gf(i). 


R(  0)  =  0,  R(  1)  =  0.  (28) 

The  above  system  of  equations  was  solved  numerically  using  PITCON  on 
the  CRAY  Y-MP/48. 

The  equations  exhibit  interesting  multiple  solutions  and  a  detailed 
discussion  of  these  solutions  and  plots  for  the  velocities  can  be  founded  in  the 
attached  paper  titled  "Multiplicity  of  solutions  in  von  Karman  flows  of  viscoelastic 
fluids". 

(ii)  Flow  of  non-Newtonian  fluids  of  the  integral  type  between  parallel  plates 
rotating  about  distinct  axes 

While  rate  type  and  differential  type  models  are  useful  in  describing  dilute 
polymeric  solutions  and  materials  with  fading  memory,  when  it  comes  to 
describing  the  behavior  of  polymeric  material  which  have  finite  memory  we  have 
to  use  integral  representations  for  the  stress  which  incorporates  the  history  of  the 
deformation  in  the  model.  The  K-BKZ  model  (cf.  Kaye  [4],  Bernstein,  Kearsley  & 
Zapas  [5])  has  proved  to  be  quite  successful  in  describing  behavior  of  a  wide  class 
of  polymeric  fluids  and  in  addition  the  model  has  been  shown  to  have  a  firm 
footing  from  the  point  of  view  of  statistical  theories  for  modeling. 

Here,  we  shall  describe  the  flow  of  a  K-BKZ  fluid  in  an  orthogonal 
rheometer  (cf.  Maxwell  &  Chartoff  [6],  Rajagopal  [7]). 

A  detailed  discussion  of  the  above  equations  for  a  special  subclass  of  K-BKZ 
fluid,  the  Wagner  fluid,  can  be  found  in  the  attached  paper  titled  "Flow  of  K-BKZ 
fluids  between  parallel  plates  rotating  about  distinct  axes:  shear-thinning  and 
inertial  effects"  while  the  flow  of  another  popular  subclass,  the  Curie  fluid  can  be 


found  in  the  attached  paper  titled  "Flow  of  viscoelastic  fluids  between  rotating 
plates  about  distinct  axes". 


The  Cauchy  stress  T  in  the  K-BKZ  fluid  has  the  structure 


T  =  -pi  +  2j'_  {ulc;‘(t)-u2c,(T)}dt. 

where 


C,(r)  =  F^(t)F,(r).  (30) 


In  (29)  U  denotes  the  strain  energy  function  for  th  e 

a  function  of  the  principal  invariants  of  Ct  (x)  and  C^1  (x): 

viscoelastic  fluid  and  is 

U  =  (/(/,, /2,r-r), 

(31) 

Ix  =  trC;1  T,  I2  =  trC,t, 

and 

U,  =  ~,  i  =  l,2. 

(32) 

ar, 

(33) 

Let  x  =  (x,  y,  z)  denote  the  position  occupied  by  the  same  particle  X  at  time  t. 

It  follows  that 

i  =-0(77  -*($)), 

(34) 

r)  =  -n(^  -f(g)). 

(36) 

?=o, 

with 

4(0  =  *.  n(0  =  )'»  and  £(*)  =  *• 

(36) 

Rajagopal  [7]  has  shown  that  the  motion  (34M36)  is  a  motion  with  constant 
principal  relative  stretch  history.  In  such  motion,  the  stress  is  determined  by  the 
first  three  Rivlin-Ericksen  tensors  Alt  A2,  and  A3. 

However,  for  the  motion  under  consideration 


Thus,  the  stress  is  given  by 


(37) 


r  =  -pl  +  /(4,,42). 

The  balance  of  linear  momentum  has  the  form: 


(38) 


if 

The  specific  constitutive  equations  determine  the  functions  hlt  h2,  and  h3. 


(39) 

(40) 

(41) 
This 


can  then  be  substituted  into  (39M41)  and  the  appropriate  partial  differential 
equations  analyzed.  Notice  that  (39M41)  are  of  second  order  and  hence  the  no-slip 
boundary  conditions  are  sufficient  for  determinacy.  The  appropriate  boundary 
conditions  for  the  velocity  field  are  (cf.  Figure  2) 


and 


u  — - Py,  v  =  fit,  w  =  0  at  z  =  0, 

2 


(43) 


u  — » v  — »  ±0°,  as  Jt,y-*±°o-  (44) 

It  follows  from  (42),  (43),  and  (34M36)  that 


/(A)  =  /(0)  =  0, 

g(h)  =  ^,  g( 0)  =  -^- 


(45) 

(46) 


In  eliminating  the  pressure  field  we  have  raised  the  order  of  the  equations.  Thus 
the  boundary  conditions  (45)  and  (46)  are  not  sufficient  to  determine  the  solution  to 
the  system  (39)  and  (41).  We  augment  the  number  of  boundary  conditions  by 
recognizing  that  the  locus  of  the  centers  of  rotation  cuts  the  plane 
z  =  0  at  some  point,  say  (/lf  t2).  Thus 


<fh'  (47) 

However,  if  we  restrict  ourselves  to  solutions  which  have  midplane  symmetry, 
then 


'(fb  <i)=°-  (48, 

For  the  rest  of  this  section  we  restrict  ourselves  to  a  discussion  of  solutions  which 
possess  midplane  symmetry.  For  the  motion  under  consideration,  a  lengthy  but 
straightforward  computation  yields 


C,(r)  =  l-^+ 


(49) 


and 


C‘(r)  =  H~[l  +  2(1  -  cir + SJ)]A -^[1  +  2(1- cir +g!)Jx!+iFA: 


+^i+£ifW2m,), 


where 


s  =  sin£l(t  ~  r),  c  *  cosQ(t- r) 


Also,  notice  that 


(50) 


(51) 


/,(/,  t  )  =  I2(t,  t)  =  3  +  2(1  -  c)(/'2  +  g'2)  s  I(Cl(t  -  t),z)  . 
It  follows  from  (39M41),  (40),  and  (50)  that 


(52) 


where 


and 


(53) 

£{-/'A(*r)+8'B(ic)}  =  pa2g, 

(54) 

<c=(r2+*T 

(55) 

A(k)  =  2 £/[3  +  2(1  -c<wl2a)ir2,a]jiflQada, 

(56) 

B(k) = 2^  {/[3+ 2(1  - cosOa)K’,a]  (1  - cos(la)da, 

(57) 

0{I,a)*U,{l,l,a)+Uj[l,l,a). 

(58) 

Lettx  and  ty  denote  the  x  and  y  components  of  the  traction  on  the  plates.  It 


follows  that 


/,(*)  =  +B(«r,  £!)/'(*)+ A(K.O)g’(h),  (59) 

»,(/i)  =  -A(ir,Q)/'(A)+A(ir,Q)*'(A)  .  (60) 

Thus  the  material  parameters  A(k,  Q)  and  B(k,  Q)  can  be  expressed  in 
terms  of  t*  and  ty  as  (cf.  Bower,  Wineman  &  Rajagopal  [8]): 


-pi*/' +'*'].  (61) 

B(K,C£)  =  -^tJ  +  tyg  ]  •  ^02) 

Dai,  Rajagopal  and  Szeri  [9]  studied  the  flow  of  a  Currie  fluid  for  which  the  strain 
energy  U  is  given  by 


U  =  t/(/„/2,s)  =  -G(s)[5/n(/  - 1)  -  9.73],  (63) 

where 

7  =  /1  +  2(/2+3.25)V2  (64) 

and 

C(s)  =  -C«-**.  (65) 

This  model  is  a  very  popular  model  for  polymeric  materials. 

In  a  recent  study  Zhang  and  Goddard  [10]  studied  the  flow  of  a  Currie  fluid 
in  an  orthogonal  rheometer  and  found  problems  with  convergence  when  values  of 
certain  parameters  were  large.  However,  Dai,  Rajagopal  and  Szeri  [9]  were  able 
to  go  well  past  the  of  the  parameters  for  which  Zhang  and  Goddard  [10]  had 
difficulty,  using  an  analytic  continuation  technique.  The  details  of  the  results  can 
be  found  in  the  paper  of  Dai,  Rajagopal  and  Szeri  [9]  appended  to  this  report.  Of 
particular  interest  is  the  fact  that  Dai,  Rajagopal  and  Szeri  [9]  were  able  to  study 


the  problem  for  Reynolds  numbers  as  high  as  10,000  and  found  the  presence  of 
interesting  boundary  layers. 

(iii)  Flow  of  non-Newtonian  fluids  due  to  torsional  and  longitudinal  oscillations 

of  a  cylinder 

There  are  many  important  applications  like  the  drilling  for  oil  where  we 
encounter  the  flow  induced  due  to  the  torsional  and  longitudinal  oscillations  of  a 
cylinder  in  a  non-Newtonian  fluid.  Here,  we  discuss  the  flow  in  the  case  of  a 
general  incompressible  simple  fluid  (cf.  Truesdell  &  Noll  [11]).  The  Cauchy  stress 
T  in  a  homogeneous  incompressible  simple  fluid  is  given  by 


T  =  -pl+MjFXt-s)]. 


(66) 


where  j  is  a  general  functional  and  Ft  (t  -  s)  is  the  relative  deformation  gradient. 

s=0 

Coleman  and  Noll  [12]  developed  a  method  for  approximating  (66)  within  the 
context  of  retarded  motions.  They  showed  that  up  to  order  four 


r  =  -pi+X  s : 

i»l 

51  = 

52  =  #i4i  +  0^4  » 

^3  =  Pi  4}  +  Pi  (A  4  +  44)  +  P^p  42)  4* 
s* =  7i4+  r2(44 + 44) + ^4 + 'J'Mi  +4  4] 

+y5(trA?)A2+  y6(trA?)Af+  y^trA^A,  +  7, (^4)4,  (71) 

where  Ax  is  as  defined  earlier  and  (cf.  Rivlin  &  Ericksen  [13]): 


(67) 

(68) 

(69) 

(70) 


4,  =  —  A.-1  +  4.1  ( grad  v)  +  ( grad  v)r4  . 


(72) 
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Using  the  above  approximation  for  the  stress,  we  studied  the  longitudinal 
and  torsional  oscillations  of  a  solid  cylinder  in  a  simple  fluid.  We  assume  a 
velocity  field  of  the  form  (cf.  Rajagopal,  Kaloni  &  Tao  [14]) 

v,  =  0,  ve  =  v(r,/),  v,  =  w{r,t)  (77) 

where  vr,  ve  and  vz  denote  the  components  of  the  velocity  in  the  r,  6  and  z 

directions.  We  assume  that  the  velocity  components  v  and  w  and  the  pressure  p 
can  be  expanded  as  a  power  series  in  Q,  for  sufficiently  small  Cl  (Cl  being  the 
frequency  of  torsional  oscillations): 


v  =  ^Q"v(,,)  +  0(£2*+1), 

*=i 

*=1 

P  =  X«V")+0(Q-+1). 


(78) 

(79) 

(80) 


Substituting  the  above  expansions  into  the  balance  of  linear  momentum  and 


equating  the  coefficients  of  Qn,  n  =  1,  ...  leads  to  a  hierarchy  of  equations  which 

can  be  solved  successively.  Here,  we  shall  not  document  these  equations  as  they 
are  exceedingly  cumbersome.  We  refer  the  reader  to  the  attached  paper  titled 
"Longitudinal  and  torsional  oscillations  of  a  solid  cylinder  in  a  simple  fluid"  for 
details  of  the  calculations. 

The  appropriate  boundary  conditions  are 

v  =  ClvcosCltee  +  Clwcos  ACl/e,,  (81) 

since  we  are  assuming  that  the  frequency  of  oscillations  in  the  circumferential 
and  longitudinal  directions  are  not  necessarily  the  same. 


Such  a  perturbation  approach  is  very  popular  and  has  been  used  by  many 
in  the  field  of  non-Newtonian  fluid  mechanics.  However,  such  an  approach  is  not 
without  drawbacks  and  these  are  also  discussed  in  some  detail  in  the  above 
mentioned  paper. 

The  study  delineates  the  effects  of  the  various  material  parameters  on  the 
wall  shear  stress.  In  general,  in  keeping  with  expectations,  an  increase  in  the 
higher  order  viscosities  like  P3  leads  to  an  increase  in  the  wall  shear  stress,  as 

also  does  an  increase  in  ft.  The  effect  of  variation  in  the  normal  stress  moduli  a! 

is  quite  dramatic.  Details  regarding  the  variations  of  the  shear  stress  T^,  v  and  w 

with  the  various  material  parameters  are  provided  in  Figures  1-12  of  the  paper 
cited  above. 

(iv)  Flow  of  non-Newtonian  fluids  in  pipes  of  varying  cross-sections 

The  flow  of  non-Newtonian  fluids  through  axially  symmetric  pipes  of 
varying  cross-sections  has  relevance  to  several  technologically  significant 
problems  in  biofluid  dynamics  and  extrusion  of  polymeric  materials.  This 
problem  has  been  studied  within  the  context  of  the  classical  linearly  viscous  fluid 
model,  and  also  specific  non-Newtonian  fluid  models  (mostly  power  law  models) 
by  several  authors. 

We  studied  the  flow  of  a  four  constant  Oldroyd  fluid  (cf.  Oldroyd  [15])  in  an 
axi-8ymmetric  pipe  using  a  perturbation  approach,  using  two  parameters:  £ 
which  describes  the  departure  of  the  cross-section  of  the  pipe  from  circularity  and 
W  (Weissenberg  number)  which  is  the  ratio  of  the  relaxation  time  to  a 
characteristic  time  scale  for  the  problem. 


The  pressure  drop  along  the  axis  of  the  tube  and  the  shear  stress  are 
calculated  up  to  second  order.  The  manner  in  which  the  elasticity  of  the  fluid 
alters  the  flow  pattern,  the  effect  of  shear-thinning  on  the  flow,  and  the  formation 
of  eddies  and  the  manner  in  which  the  flow  separates  are  all  studied  in  detail. 
When  appropriate  parameters  are  set  to  zero,  the  results  reduce  to  the  results 
established  for  the  Navier-Stokes  fluid.  We  find  that  elastic  effects  can  initiate 
separation  at  much  lower  Reynolds  number  than  for  a  Newtonian  fluid. 
Moreover,  the  wall  shear  stress  can  change  dramatically  with  changes  in  some  of 
the  non-Newtonian  parameters.  For  instance,  increasing  the  retardation 
parameter  decreases  the  wall  shear  stress,  while  increasing  another  non- 
Newtonian  parameter  (a)  which  appears  in  the  model  can  significantly  increase 
the  wall  shear  stress. 

The  Cauchy  stress  T  in  a  four  constant  Oldroyd  fluid  is  related  to  the  fluid 
motion  in  the  following  manner  (cf.  Oldroyd  [  ]): 


T  +  XA 


(82) 


where  D/Dt  represents  the  co-rotational  or  Jaumann  time  derivative,  given  by  (for 
any  symmetric  second  order  tensor  B) 


and 


DB 

DT 


dB 

dt 


=  [grad  B\u  +  WB  -  BW, 


(83) 


A  =  Vu  =  (Vu)r,  2W  =  Vu-(Vu)r. 


(84) 


Also,  t|q  is  the  viscosity  of  the  fluid  at  the  zero  shear  rate,  Xj,  is  the  relaxation 


time,  X2  is  the  retardation  time  and  a  is  a  parameter  whose  value  lies  between  0 

and  1.  When  a  =  0,  the  above  model  is  equivalent  to  Jeffery’s  version  of  the 

Oldroyd  model  whereas  when  a  =  1  it  reduces  to  the  limiting  form  of  Walter’s  four 
constant  version  of  Oldroyd's  model.  In  viscoelastic  materials  it  seems 
reasonable  to  assume  that  0  <  a  <  1.  In  this  range  predicts  shear  rate  dependent 
viscosity,  normal  stress  effects  in  shear  flow  and  a  reasonable  behavior  in 
elongatial  flow. 

We  assume  the  flow  to  be  axisymmetric  so  that  the  components  of  the 
velocity  (u,  v,  w)  in  (x,  r,  0)  direction  satisfy  the  constraint  of  incompressibility: 

du  dv  v  _ 

— +  —  +  -  =  0. 

obc  dr  r  (86) 

We  introduce  the  axi-symmetric  stream  function  \j/  through 


¥ 

u  =  —  and  v  =  — 
r  r 

and  the  vorticity  component  Q  through 


(86) 


Q  =  ur-vx  = 


(87) 


where  the  subscript  denotes  partial  differention  with  respect  to  that  variable. 
The  appropriate  boundary  condition  on  3D  are 


on 


r  =  a(x) 


V  =  0,  v  =  0,  and  ur  =  0  on  r  =  0 . 


(88) 


Equations  (88)1,2  are  a  consequence  of  the  adherence  boundary  condition  and 

imply  that  there  is  no  fluid  motion  either  tangential  or  normal  to  the  wall.  The 
conditions  (88)3,4  are  essentially  symmetry  conditions  on  the  axis  of  the  tube.  The 

constant  flow  rate  taking  place  through  any  cross-section  of  the  tube  is  given  by 


ff 

Jo  Jo 


rudddr  =  2  JtQ  =  constant . 


The  boundary  conditions  (88)1.4  can  be  expressed  in  terms  of  by  entering  (88) 
into  (86): 


4^  =  0,  '¥  =  '¥0  on  r  =  a(x) 
4?  =  0  on  r  =  0 


PI 


-»  0  as  r  — » 0 . 


We  now  proceed  to  non-dimensionalize  the  governing  equations.  Let 


u  =  ^u,  v  =  ^fv,  f  =  A  =  ^fA,  W  = 

Oo  <*o  %  c$  £ 

P  =  V  =  n  =  ('P0ja0>(z,x*,£), 

a(x,€)  =  a0j(a|^o),  (0<e«l),  x* ^exfa,  z  =  r|ao.  (91) 

where  p  is  the  pressure,  and  ao  is  a  constant  characteristic  radius  of  the  tube. 
The  function  s(ex|ao)  is  such  that  in  the  limit  e  ->  0,  the  tube  is  of  constant  radius, 

and  the  variation  of  a  with  respect  to  the  axial  coordinate  x  depends  upon  ex 
rather  than  x  alone.  Here  z  is  the  normalized  radial  coordinate  and  x*  is  a 


"slowly  varying"  normalized  axial  coordinate.  On  substituting  (91)  into  (82)  we 
obtain  the  dimensionless  form  of  the  constitutive  equation  as 


where 


T  =  2A  +  2ex  =  WA  -  AW  -  2a(42)  j 

-Wi{^  +  WT-TW~a(AT+TA)}> 


Wi  =  — p3-  ( Weissenberg  number) 


(92) 


and 


1  v 

The  dimensionless  form  of  the  equation  of  motion  in  the  absence  of  body 
forces  is  given  by 


dv 

div  T  -  grad  p  =  Re — , 
dt 


(93) 


where  Re  =  ^(/flo  v  is  a  Reynolds  number. 

On  substituting  (91)  into  (90)  we  get  the  dimensionless  boundary  conditions 
as 


0=0  on  :  =0,  <p,./z->0  as  z->  0,  {4Jz)t-*0  as  z  ->0 .  (94) 

We  solve  the  system  of  equations  (92)-(94)  for  small  Weissenberg  and  Reynolds 
number.  The  previous  low  Reynolds  number  solutions  for  the  Navier-Stokes 


equations  can  be  readily  retrieved  from  our  solution  and  our  results  reduce  to  the 
results  of  Manton  [16]  when  the  Weissenberg  number  Wi  =  0. 

First,  assuming  that  Re  is  fixed  and  is  of  order  unity,  a  solution  is  obtained 
in  the  form  of  a  power  series  in  Wi  for  the  quantities  <{>,  to,  u,  v,  q,  T,  e,  and  w  (cf. 
Kasivishvanathan,  Kaloni  &  Rajagopal  [17]): 

A  =  A0  +  WiAi  +  (Wi)2A2+...  (95) 

where  A  in  turn  stands  for  <j>,  co,  u,  v,  q,  T,  e  and  w. 

Details  of  the  results  can  be  found  in  the  enclosed  paper  titled  "Flow  of  a 
non-Newtonian  fluid  through  axi-symmetric  pipes  of  varying  cross-sections". 
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